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We consider the lightly doped cuprates Yo.grCao.osBaCuOe.os and La 2 -ccSra;Cu 04 (with x — 

0.02,0.04), where the presence of a fluctuating nematic state has often been proposed as a precur¬ 
sor of the stripe (or, more generically, charge-density wave) phase, which sets in at higher doping. 

We phenomenologically assume a quantum critical character for the longitudinal and transverse 
nematic, and for the charge-ordering fluctuations, and investigate the effects of these fluctuations 
in Raman spectra. We find that the longitudinal nematic fluctuations peaked at zero transferred 
momentum account well for the anomalous Raman absorption observed in these systems in the B 2 g 
channel, while the absence of such effect in the Big channel may be due to the overall suppression of 
Raman response at low frequencies, associated with the pseudogap. While in Yo.grCao.osBaCuOe.os 
the low-frequency lineshape is fully accounted by longitudinal nematic collective modes alone, in 
La 2 -ccSrxCu 04 also charge-ordering modes with finite characteristic wavevector are needed to re¬ 
produce the shoulders observed in the Raman response. This different involvement of the nearly 
critical modes in the two materials suggests a different evolution of the nematic state at very low 
doping into the nearly charge-ordered state at higher doping. 

PACS numbers: 74.72.-h, 74.25.nd, 75.25.Dk, 74.40.Kb 


I. INTRODUCTION 

Growing experimental and theoretical evidence indi¬ 
cates that (stripe-like) charge ordering (CO) [IH3], pos¬ 
sibly related to a hidden charge-density-wave quantum 
critical point near optimal doping mu, plays a role 
in determining the unconventional properties of super¬ 
conducting cuprates. Charge ordered textures were as¬ 
sessed by neutron scattering experiments in La cuprates, 
codoped with Nd [SUIT]. Ba [12], or Eu [T3|, and con¬ 
firmed also by soft resonant x-ray scattering mug. The 
occurrence of stripe-like charge- and spin-density waves 
in other cuprates is supported by the similarities of the 
non-codoped and codoped La cuprates in the spin chan¬ 
nel, e.g., the doping dependence of the low-energy incom¬ 
mensurability m, and the high-energy magnon spectra 
in La 2 -a;Baa,Cu 04 [TT] . La2-ccSra,Cu04 (LSCO) [T^ . and 
YBaCuOe+p |19I2Q| . These features are well described in 
terms of striped ground states EBEg. CO in cuprates, 
possibly with fluctuating character, was also confirmed 
by EX AES [24], NMR experiments [25126) . scanning tun¬ 
neling spectroscopy [27H29] , and resonant x-ray measure¬ 
ments [301433] . A recent theoretical analysis of Raman 
spectra in LSCO [34] showed that nearly critical spin 
and charge fluctuations coexist at intermediate and high 
doping. This coexistence also accounts [35] for the spe¬ 
cific momentum, energy and doping dependence of the 
single-particle anomalies, the so-called kinks and water¬ 
falls, observed in photoemission spectra [36] . 

The above facts, support the occurrence and relevance 
of (fluctuating) stripes in cuprates and raise the ques¬ 
tion about their precursors at very low doping [37]. The 
experimental evidence of rotational symmetry breaking 


[201 [MSI] points towards nematic order, although it is 
not yet clear whether this order arises from a melted 
stripe state [42], from incipient unidirectional fluctuating 
stripes [33], or from an unrelated (i-wave-type nematic or¬ 
der which preserves translational symmetry [44]. On the 
theoretical side, it was recently proposed that a ferrone- 
matic state occurs at very low doping, formed by stripe 
segments without positional order m- These segments 
are oriented because they sustain a vortex and an antivor¬ 
tex of the antiferromagnetic order at their extremes, and 
break rotational and inversion symmetry. This phase has 
no order in the charge sector, but induces incommensu¬ 
rate peaks in excellent agreement with experiments in 
LSCO [46]. Recent Monte Carlo calculations m showed 
that, lowering the temperature, the ferronematic state 
turns into a ferrosmectic state, where the segments have 
a typical lateral distance ic^ corresponding to CO with 
a characteristic wavevector (with jqd ^ 1/^c)- The 
segments thus appear as the natural precursors of stripes. 

It is therefore important to assess nematic order in 
cuprates. The aim of the present work is to identify 
the signatures of nematic fluctuations in Raman scat¬ 
tering. This is a bulk (nearly surface-insensitive) probe 
and measures a response function analogous to that of 
optical conductivity [48]. However, while the latter av¬ 
erages over the Brillouin zone (BZ), different polariza¬ 
tions of the incoming and outgoing photons weight dif¬ 
ferent parts of the BZ in Raman scattering m , intro¬ 
ducing specific form factors. It turns out that the so- 
called Big and B 2 g channels are the most relevant to 
extract the contributions of collective modes (CMs) in 
cuprates. We already investigated how these form fac¬ 
tors can be exploited to identify the contributions of dif- 
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ferent (e.g., charge and spin) critical CMs, based on their 
different finite wavevectors |34I5QI53H55] . There are two 
classes of CM contributions. In one class, the CMs dress 
the fermion quasiparticles, introducing self-energy and 
vertex corrections, which affect the Raman spectra up 
to substantial fractions of eV |34I53I54| . In the other 
class, the excitation of pairs of CMs m, reminiscent of 
the Aslamazov-Larkin (AL) paraconductive fluctuations 
near the metal-superconductor transition (see Fig. §,af- 
fects mainly the low frequency part of the spectrum and 
produces an anomalous absorption up to few hundreds 
of cm“^, as indeed observed, e.g., in LSCO [56]. The 
analysis for LSCO [50| was based on CO CMs with fi¬ 
nite wavevector qc, while the role of spin CMs was ruled 
out by symmetry arguments. At moderate doping, the 
value qc ~ (±7r/2, 0), (0, ±7r/2) was deduced from inelas¬ 
tic neutron scattering as the double of the wavevector of 
spin incommensuration m, within the stripe scheme (we 
use hereafter a square unit cell on the Cu 02 planes, with 
lattice spacing a = 1). By symmetry arguments, and in 
agreement with experiments, fluctuations with such qc 
give rise to an anomalous absorption in the Big channel 
only. A rotated qc ~ 27r(±2x, ±2x) occurs for x < 0.05, 
|46| . making the anomalous Raman absorption show in 
the B 2 g channel only, consistent with the experiments. 
However, a similar anomalous absorption in the B 2 g 
channel is observed in Yi_^CayBa2Cu306+cc (YBCO) for 
doping p{x^y) between 0.01 and 0.06 (52). Recent mea¬ 
surements [38] do not support the rotation of the spin 
modulation vector in YBCO, at least down to p = 0.05, 
and the extrapolation of the available data indicates that 
spin incommensuration disappears for p ^ 0.02 — 0.03, 
while CO seems to disappear for p < 0.08 [33]. Thus, if 
only CO fluctuations were to play a role, the anomalous 
peak observed in YBCO in the B 2 g channel would be un¬ 
explained. Furthermore, CO CMs yield in LSCO spectra 
that are fully satisfactory at x = 0.1, but less convincing 
at X = 0.02, where the experimental lineshape seems to 
have a composite character, with a main peak accompa¬ 
nied by a shoulder at slightly higher frequencies. This 
suggests the presence of two CMs contributing to the 
anomalous absorption in the B 2 g channel at low doping 
in LSCO and raises the question about the nature of the 
additional CM. The uncertain situation with YBCO and 
the compositeness of the LSCO spectra call for a critical 
revision of the results of Ref. m- 

The above mentioned evidences for nematic order make 
it natural to inquire whether the anomalous Raman ab¬ 
sorption observed in underdoped cuprates might be due 
to nematic fluctuations (not considered in Ref. [50]), pos¬ 
sibly mixed with CO fluctuations (in LSCO). Therefore, 
within the same formal scheme of Ref. m, we include 
here the contribution of nematic fluctuations. We find 
indeed that at low doping the observed anomalous ab¬ 
sorption can be due to the excitation of long-wavelength 
overdamped nematic fluctuations with longitudinal char¬ 
acter, whose strong dynamics is apt to reproduce the ob¬ 
served lineshape. While in strongly underdoped YBCO 


this is enough, in LSCO a secondary CM with finite char¬ 
acteristic wavevector, which we identify with the CO CM, 
is needed to better represent the lineshape. The doping 
dependence of the lineshape in LSCO indicates that there 
is an evolution from a dominating NCM towards a major 
relevance of the CO CM, upon increasing doping. 

The scheme of the paper is the following: In Sec.|TI|we 
introduce a phenomenological model of fermion quasipar¬ 
ticles coupled to nearly critical CO CMs and NCMs in 
underdoped cuprates. Then, we proceed with the theo¬ 
retical calculation of the Raman response due to these 
CMs (Secs.l^l^ and lira. In Sec.lm, we compare 


the theoretical results with available Raman spectra for 
underdoped YBCO and LSCO. Sec.|^ contains our fi¬ 
nal remarks and conclusions. Appendix A contains some 
details about the calculations of the Feynman diagrams 
involved in the anomalous Raman response. Details of 
the fitting procedure are found in Appendix B, while a 
discussion on the role of the pseudogap in the fermionic 
spectrum is found in Appendix C. 


II. THE FERMION-COLLECTIVE MODE 
MODEL AND THE RAMAN RESPONSE 

A. The fermion-collective mode model 

We consider a phenomenological model where, simi¬ 
larly to the electron-phonon coupling, electrons are cou¬ 
pled to NCMs or CO CMs. This approach relies on 
the presence of fermion quasiparticles. This assumption, 
which is natural in the metallic phase of cuprates, is still 
justified in the strongly underdoped phase, where angle 
resolved photoemission [S71 and [55] and transport exper¬ 
iments [59] and EO] highlight the presence of fermionic 
low-energy states (the so-called Fermi arcs) with a sub¬ 
stantial mobility, indicating that fermion quasiparticles 
still survive in this “difficult habitat”. Thus, we adopt 
the Hamiltonian 

^ = Y (1) 

k,cr k;,q,cr A 

where (cka) creates (annihilates) a fermion quasipar¬ 
ticle with momentum k and spin projection a, and 
is the fermion dispersion on the Cu02 planes of LSCO 
or YBCO (measured with respect to the chemical po¬ 
tential). Its specific form is rather immaterial for our 
analysis, once the generic shape of the Fermi surface of 
cuprates is taken into account. The index A labels trans¬ 
verse (A = t) or longitudinal (A = i) nematic fluctu¬ 
ations Ensi], and charge fluctuations (A = c), repre¬ 
sented by the boson fields The quasiparticles couple 
to NCMs via g\{k,q) = with ^ = cos(2(^k,q) 

and q = sin( 2 (pk,q )5 where (pk,q is the angle between 
k and q (see, e.g., Eqs. (2.3) in Ref. [62]). The CO CM 
has instead a finite characteristic wavevector qc and cou¬ 
ples to the fermion quasiparticle via a weakly momentum 
dependent coupling ^c(k,q) « Pc (i.e., ^ ^ 1). 
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We assume that these CMs are near an instability and 
their propagators take the standard Gaussian form, valid 
within a Landau-Wilson approach, and already adopted 
for models of fermion quasiparticles coupled to nearly 
critical charge [4] and spin CMs in cuprates. 

As customary in quantum critical phenomena, different 
damping processes may lead to different dynamical criti¬ 
cal exponents 2 ;, relating the divergent correlation length 
^ and time scale r oc In the case of the nematic 
instability a multiscale criticality occurs due to the dif¬ 
ferent dynamics of transverse and longitudinal fluctua¬ 
tions IK]. The longitudinal fluctuations are Landau- 
overdamped, decay in particle-hole pairs acquiring a dy¬ 
namical exponent Z£ = 3, and their propagator is 

I 12 I I I /I I I 2 /O ’ 

where ujn is a boson Matsubara frequency and wavevec- 
tors q are henceforth assumed dimensionless and mea¬ 
sured in units of inverse lattice spacing a~^ (when 
needed, conventional units are restored in our formulas 
by replacing q with aq). Apart from the term oc , this 
propagator is the same as that in Eq. (2.14) of Ref. [62l 
Transverse fluctuations have instead = 2, and their 
propagator is [see, e.g., Eq. (2.15) in Ref. [62] 


Both propagators, in the static limit {ujn = 0), are peaked 
at q = 0. Similarly, the nearly critical CO CM has a 
dynamical critical index Zc = 2, with propagator [see, 
e.g., Eq. (1) in Ref. [51] or Eq. (2) in Ref. |3i] 


T>c(q, ^n) 


1 

'^c T Cc|q qcp T |^n| T ^n/^c 


(4) 


peaked at a finite wavevector qc (actually, at the whole 
star of equivalent wavevectors). This circumstance al¬ 
lows to reabsorb a factor \q_c\‘^ in the definition of the 
parameter and marks the difference with respect to 
the propagator of the transverse NCMs, Eq. In the 
doping regime we are considering, qc is directed along 
the diagonals of the BZ in LSCO with x < 0.05 [16]. 
According to the discussion in Sec.[T] we consider instead 
that CO is absent in YBCO with p ^ 0.015. 

In Eqs. ([2][^, the parameters cx set the curvature at 
the bottom of the CM dispersions, whereas the param¬ 
eters Clx set high-frequency cutoffs. The low-frequency 
scales mx are proportional to the inverse squared correla¬ 
tion lengths thus being the relevant parameters that 
measure the distance from criticality. 


B. The fermionic loop in the Raman response 

Our theoretical analysis is based on the calculation of 
the Raman response represented by the Eeynman dia¬ 
grams of Eig.[^ (more details are given in Appendix [A|). 



Figure 1. Diagrammatic representation of the Raman re¬ 
sponse due to the excitation of two CMs. The grey dots 
represent the 'jBig or 'yB 2 g form factors. The solid lines are 
the propagators of the fermion quasiparticles in the fermionic 
loops, the w avy lin es represent the NCM or CO CM propa¬ 
gators [Eq.s which are coupled to the quasiparticles 

by the coupling functions ^A(k, q) (solid dots). 


The first step is to calculate the sum of the fermionic 
loops with attached direct and crossed boson lines (see 
top and bottom diagrams in Eig.[^ 

Ay(q,z/,,w„) = ry]7j(k)5rA(k,q)fl,,(k,-q) 

k,n 

x[G(k + q, ^rn ) + G(k + q, + vi)] 

xG(k,en)G(k,en + i^i), (5) 

where T is the temperature, i = Big,B 2 g labels the 
form factors, 7Bi^(k) = cos{ky) — cos{kx) and 7^2^ (k) = 
sm{kx) sm{ky) [S], is the external Matsubara fre¬ 
quency which, once analytically continued, represents the 
frequency shift between the incoming and the scattered 
photons, ujm is the Matsubara frequency of one of the bo¬ 
son propagators in Eig.[^(the other carries cj^ + z//), is 
the fermion frequency to be summed over in the fermionic 
loop, and G(k, e^) = (ie^ — ^k)~^ is the fermion quasi¬ 
particle propagator. In Eq. ^ we exploited the parity of 
G(k, Cn), 7i(k), and 5 fA(k, q)^^(k, -q) with respect to k. 

The dependence of the loop on the CM indexes A and 
T] is diagonal: the CO CM cannot mix with the NCMs, 
having a finite characteristic wavevector, and the i and t 
NCMs cannot mix, because the product of 5 f^(k, q) and 
^t(k, —q), each depending only on the angle between k 
and q and having a different parity, averages to zero when 
summed with respect to k. This fact entails a selection 
rule stating that the two NCMs attached to the same 
fermionic loop must be either longitudinal or transverse. 
The average over the Eermi surface of two couplings with 
the same NCM yields a result that is weakly dependent 
on q and can be safely approximated to a constant that 
can be reabsorbed in the definition of the dimensional 
coupling gx. Thus 
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Summing over the fermion frequencies, one obtains the 
general expression 


= 2^ 

k 


Ti(^)^/k ^m{^m “1“ 


where A/k = /(^k+q) - /(^k), A^k = ^k+q - ^k, and 
f{z) = (e^/^ + 1)“^ is the Fermi function. 

The next steps are different in the case of NCMs (with 
characteristic wavevectors q ^ 0) and of CO CMs (with 
finite characteristic wavevectors qc), and will be dealt 
with in Sec. 1 B l| and Sec. |IIB 2} respectively. 


1. The fermionic loop for NCMs 


To proceed with the calculation of the fermionic loop 
in the case of NCMs, we consider that the main features 
of the boson propagators and (§ are their poles at 
small momenta g = |q| and even smaller frequencies, be¬ 
cause of their dynamics with zi = 3 {uj ^ q^) or Zt = 2 
{uj ^ q^). Thus, expanding the above result for small fre¬ 
quencies and keeping the lowest order in the Matsubara 
frequencies and + ui , one obtains 


k 


7i(k)A/k 


2 '^ 7i(k) df{^k) 


The summation on k can be transformed into a two- 
dimensional integral, yielding 

Ai(q) dkdeSik - (6) 

where 0 is the angle between the wavevector k and 
the X axis in reciprocal space, and M is the quasipar¬ 
ticle effective mass. By noticing that the form factor 
7 i(k), calculated on the Fermi surface, depends weakly 
on k = |k| while it substantially depends on one can 
write ^Big{k,0) ^ cos(26>) = ^Big{0) and ^B 2 g{k, 0 ) « 
sin(26>) = jB 2 g{^)‘ When expanding A^k one has to keep 
track of the inverse band curvature M (otherwise the in¬ 
tegral vanishes). The limit |q| ^ 0 can then be taken, 
and the final result is that the fermionic loop depends 
only on the angle (f between q and the x axis. This de¬ 
pendence can be made explicit observing that the denom¬ 
inator A^k in Eq. depends on the cosine of the angle 
0 — (j) between k and q. Shifting the variable 0 — (p ^ 0^ 
one is left with ^i{0 + p) in the numerator. Expanding, 
one has jBlg{0 + 0 ) = TBlg{0)jBlg{<P) — jB2g{0)^B2g{(p) 

andjB2g{0 + (t)) = ^B2g{0)jBlg{(t>)^lBlg{0)^B2g{(t>)^ The 
integral with respect to 0 of the terms with ^B 2 g{ 0 ) van¬ 
ishes by symmetry. Thus, we finally obtain 

~ (</>)• (7) 

This result, which is crucial in our development, implies 
that the original form factor ^i{0) coupling the fermion 


quasiparticles to the incoming and outgoing photons in 
the Raman vertex, in the integrated form of the loops, 
is translated into a direct coupling of the photons to the 
NCMs with the same form factor 7i(0). 


2. The fermionic loop for the CO CMs 


The fermionic loop for the CO CMs has been calculated 
in Ref. [50], and we recall here the main results. The main 
difference with respect to the calculation of Sec. |IIB fj is 
that the propagator is peaked at finite wavevectors qc- 
Then, the sum over k in Eq. ^ is now dominated by the 
neighborhood of the points along the Eermi surface where 
^k = ^k+qc 5 i -^-5 so-called hot spots. The results is 


Ai(qc)«^ln 




W- 


^ li,HS 

v%g sin aHS 


(8) 


where W± are the upper and lower cutoffs for the lin¬ 
earized band dispersion at the hot spot, while ji^HS aiid 
vhs are, respectively, the Raman form factor and the 
Eermi velocity at k = ki:/^, and ans is the angle be¬ 
tween the Eermi velocities at the two hot spots connected 
by the given qc. Eor qc along high symmetry directions 
(i.e., the axes and the diagonals) of the BZ, the moduli 
of the Eermi velocities in ]<.hs and ]<.hs + Qc are equal. 
As pointed out in Ref. m, summing over k at fixed qc, 
various different hot spots are visited, where, due to the 
above-mentioned symmetry, the form factors can have 
pairwise equal magnitude and equal or opposite signs. 
As a consequence, the terms in the above hot-spot sum¬ 
mation can add or cancel each other. This induces a “se¬ 
lection rule” which, in the case pertinent to the strongly 
underdoped LSCO, where qc is short and directed along 
the (±1,±1) directions, leads to finite B 2 g vertex loops, 
while the Big vertex loops vanish by symmetry. 


C. The Raman response 

Eew considerations are now in order. Eirst of all, the 
NCM propagators, Eqs. and do not depend on 
the angle p and therefore the product of the two fermionic 
loops entering the diagrams of Eig.[2 only introduces a 
multiplicative constant factor, which can be enclosed in 
the overall intensity of the Raman response. However, we 
emphasize that the p integration, to be performed when 
the summation over q is carried out, introduces an im¬ 
portant selection rule: The fermionic loops with attached 
Raman vertices, enter pairwise in the response diagrams 
of Eig.fBand must both be of the same symmetry. Big or 
B 2 g. Similarly, the CO propagator 0 depends only on 
the magnitude of the deviation of q from qc- In this case 
the Raman response is given by a first summation on all 
the possible qc of [A^(qc)]^ and an internal integral over 
|q — qd of two CO CM propagators. 
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Thus, both for the two nematic CMs and for the CO 
CM, the sum of the two diagrams of Fig.reads 


q 

Xi,\iyi) = Ki^xT'Y' / Aq^Dx{q,u:n)Dx{q,i^n + vi), 

„ Jo r 


where g ^ 1 is the momentum cutoff, its precise value 
being re-absorbable in a multiplicative rescaling of the 
parameters of the CM propagator, q = |q| for the 
NCMs, and q = |q — qd for tbe CO CMs. The fac¬ 
tor Ki^x comes from the product of two fermionic loops 
and is proportional to g\ (each loop carrying two 
fermion-CM coupling constants). For NCMs, Ki^x = 

with \ = e or t, and ([7i((/>)]^) 
is the angular average of the square of the function 7 ^ ( 0 ) 
that appears in Eq. Q, whereas for the CO CM we have 
[cf. Eq. 0] Ki^c = SqJAi(qc)]^, that vanishes in 
the Big {B 2 g) channel for diagonal (vertical/horizontal) 
qc- This “selection rule” is the only place where the finite 
wavevector of the CO CM plays a role within our nearly 
critical theory of Raman absorption. This selection rule 
is instead absent in the case of the NCMs, which are 
peaked at q = 0 . 

The analytic continuation to real frequencies iui 
uj ^ id and the use of the spectral representation of the 
boson propagators finally yield the Raman response 

/ +00 pi 

dz [b(z-) - b{z+)] / dg 
-00 J 0 

X 2qFx{z+,q)Fx{z-,q), (9) 

where b{z) = (e^/^ — 1)“^ is the Bose function, and we 
performed the customary symmetrization z ^ z — ^ = 
Z- and z^uj which makes explicit the fact 

that Eq. 0 is an odd function of uj. The constant multi¬ 
plicative prefactors, including those transforming the Ra¬ 
man susceptibility into the measured Raman response, 
are reabsorbed in the parameters Unfortunately 

a fully analytical expression for Ai^x cannot be given, 
the Raman response being affected by resonance effects 
that prevent even order-of-magnitude estimates. How¬ 
ever, whenever we studied the contribution of critical 
CMs in situations where the prefactors can be explicitly 
calculated, like optical conductivity mm, or angle re¬ 
solved photoemission spectra m, we always found that 
the dimensionless coupling constants are of order one, in 
a regime of moderate coupling. 

The spectral density of the longitudinal NCMs is 


Ft{z,q) 


q. 


{mi + ciq^ - 



while the spectral density of the transverse NCMs is 
Ftiz.q) = - - , s2 - • 



CO [cm 


Figure 2. Schematic representation of the AL-like Raman re¬ 
sponse of the three different CMs: longitudinal L-NCM (black 
curve), transverse (T) NCM (green curve), and CO CM (red 
curve). Inset: AL-like Raman response from different CMs, 
but taken with the same nearly critical set of parameters 
{mx = 17 cm“^, Ua = 70 cm“^, ca = 3.16) The ampli¬ 
tudes are instead rescaled by factors of order one to bring all 
responses to a common maximal height for easier comparison. 
The coloring of the lines is the same as in the main panel. 


Finally, for the CO CMs we find 

Fciz,q) = - . 

{me + ^ j + 

The anomalous peak of the Raman response, both in 
LSCO and YBCO, is strongly temperature dependent, 
it shrinks and softens upon reducing T. This behavior 
is naturally encoded in the temperature dependence of 
the mass mx of the CMs. In general, the low-frequency 
scale mx controls the slope of the Raman response, while 
the scales uui ^ xjmx^x and UJ 2 ^ \/(^a + ca)Ua set 
the frequency window over which the spectral function 
of the corresponding CM is sizable. However, the differ¬ 
ent dynamical properties and values of the parameters of 
the CMs mirror into different shapes of the AL-like Ra¬ 
man responses, which are schematically represented in 
Fig-H We point out that the curves displayed in this 
panel do not exhaust all the possible regimes of param¬ 
eters, and only represent the corresponding CM in the 
regime where, after a thorough analysis, they were found 
to better reproduce the various features of the Raman 
response. The inset of the same figure reports instead 
the behavior of Raman absorption spectra (from the AL 
processes) due to the various CMs. In this inset, while 
we rescale the height to bring all responses to the same 
maximal height, we use the same nearly critical set of pa¬ 
rameters to highlight the differences arising purely from 
the different form of the propagators and dynamical crit¬ 
ical index Apparently the shape of the spectra is quite 
similar, but the behavior upon changing the mass is dif- 
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ferent on a quantitative level. In particular we found that 
the z = 2 propagators shift the position of the maxima 
upon reducing m more rapidly than the NCM 2 : = 3 prop¬ 
agator. Si nce we apply a strict fitting protocol (see below 
in Section [lIl| ) , which fixes all parameters and follows the 
temperature evolutions of the main peaks by only chang¬ 
ing m, these different behavior affects in a substantial 
way the accuracy of the fits. An inspection of Fig. 3 in 
Ref. [50] shows that the fits with a z = 2 CO-CM are not 
very accurate at low temperatures. Instead we will see 
in the next Section that the 2 : = 3 NCM does a much 
better job within the adopted strict fitting protocol and 
therefore it will be considered henceforth as the primary 
(i.e., most critical) CM. The additional shoulder in the 
spectra of LSCO, is instead better reproduced by the CO 
curve in Fig.than by the broader T-NCM curve, when 
both CMs are taken in the regime of parameters apt to 
describe this spectral feature. Therefore at these doping 
levels the CO CM acts as the secondary CM in LSCO. 


III. RESULTS 

A. Raman absorption in Y0.97Ca0.03Ba2Cu3O6.05 

An anomalous Raman absorption at low frequencies, 
up to few hundreds of cm“^, is experimentally found in 
the B 2 g channel in lightly doped YBCO with p < 0.05 
[52]. Since, however, the whole spectra also display broad 
absorptions up to electronic energy scales, we first extract 
the specific anomalous low-frequency contributions. To 
this purpose, we subtract from the low-temperature spec¬ 
tra the spectra obtained at the highest measured tem¬ 
perature. This subtraction is delicate because at tem¬ 
peratures below about 150-200 K, the spectra are char¬ 
acterized by the formation of a pseudogap over a fre¬ 
quency range of several hundreds of cm“^, which reduces 
the electronic background. Then, the simple subtraction 
leads to regions of negative absorptions, which are ob¬ 
viously meaningless. In the Appendix B we provide the 
detailed procedure adopted to circumvent this drawback. 
In Fig.[^ the data processed according to the previous 
procedure, are shown for p ^ 0.015. 

The experimental lineshape clearly resembles the L- 
NCM spectrum in Fig.|^ which is narrow due to the 
Zi = 3 damped dynamics of the corresponding CM, 
whose temperature dependence is ruled by the mass 
m£. Indeed, the data in Fig.are best fitted with the 
only contribution of longitudinal NCMs. In the spirit of 
our nearly-critical approach, we only adjust their mass 
m^(T), while keeping all other parameters (i.e., the high- 
frequency cutoffs of the CM propagator, the q coeffi¬ 
cients, and the overall intensity coefficient A^) fixed at 
all temperatures. This strict procedure was already suc¬ 
cessfully adopted in Ref. m and seems to us the most 
suitable to pinpoint the quantum nearly-critical charac¬ 
ter of the collective excitations responsible for the anoma¬ 
lous Raman absorption. The fits with this restricted 



Figure 3. Subtracted experimental Raman absorption spectra 
in the B 2 g channel, at various temperatures, for YBCO at 
p ^ 0.015 (symbols). The theoretical fits (solid lines) consider 
the contribution of the longitudinal NCM only. The fitting 
parameters are q = 0.63 cm“^, fli = 110 cm“^, Ai = 5.0 
(a.u.). The inset reports the temperature dependence of the 
mass of the longitudinal NCM (black circles). 


procedure turn out to be quite good. Of course they 
could be further improved if this constrained procedure 
were relaxed. The fits reproduce well the lineshapes and 
the strong temperature dependence of the peak, encoded 
in the rapid decrease of the mass with temperature, as 
shown in the inset of Fig.[^ From this inset, it is ev¬ 
ident that mi{T) decreases with T. Its linear extrapo¬ 
lation starting from high temperature should vanish at 
some finite critical temperature for the onset of nematic- 
ity (~ 125 K), if static order would occur. However, at 
lower temperatures, the mass seems instead to saturate, 
likely indicating that nematic order stays short-ranged 
and dynamic. 


B. Raman absorption in La 2 -ccSrccCu 04 

Figs.j^ and report the experimental Raman spec¬ 
tra in the B 2 g channel, for LSCO samples at doping 
X = 0.02 and x = 0.04 and various temperatures. The 
raw data were again processed according to the proce¬ 
dure described in the Appendix B. As mentioned above, 
the anomalous Raman absorption observed in LSCO is 
characterized by a lineshape that is more complex than 
in YBCO, and displays a peculiar shoulder or, at low T, 
even a secondary peak, see Fig.[^ The anomalous peak 
and the shoulder (or secondary peak) both depend on 
temperature, but their frequency and intensity are not 
simply related by constant multiplicative factors. The 
shoulder (or secondary peak) becomes stronger with in¬ 
creasing doping. This indicates that the excitations re¬ 
sponsible for this absorption have a distinct dynamics. 
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Figure 4. Subtracted experimental Raman absorption spectra 
in the B 2 g channel, at various temperatures, for LSCO at 
X = 0.02 (symbols). The theoretical fits (solid lines) consider 
the contribution of the longitudinal NCM and of the CO CM, 
with Q = 3.16 cm“^, Cc = 333 cm“^, = 70 cm“^, — 8.3 

(a.u.). The other fitting parameters are reported in Fig.[^ 
The inset reports the temperature dependence of the mass of 
the longitudinal NCM (black circles) and of the CO CM (red 
squares). 


Again these absorptions are described by the AL-like 
processes (direct and crossed, see Fig.[^. Owing to 
the selection rules found in Sec.[^ the response due to 
two (or more) CMs is the sum of the responses associ¬ 
ated with each individual CM. As already mentioned, 
our thorough analysis showed that the primary anoma¬ 
lous absorption should be attributed to the longitudi¬ 
nal NCM, which has the stronger dynamical behavior. 
Within our context, the transverse NCM and the CO 
CM are the two candidates for the shoulder (or secondary 
peak). Looking at the lineshape of the two CMs reported 
in Fig.[^ it is easy to convince oneself that the best choice 
for a good fit is the CO CM, due to its much more pro¬ 
nounced peaked form at intermediate frequency. We also 
attempted a fit with the transverse NCM. At x = 0.02 
we obtained a reasonable fit taking a very large and al¬ 
most temperature independent CM mass, which is hardly 
compatible with our assumption of nearly critical CMs. 
Moreover, at x = 0.04, when the shoulder evolves into a 
secondary peak, the attempt failed completely. Thus, we 
ruled out a contribution of transverse NCMs. 

Again, having attributed the main peak to the more 
critical longitudinal NCM, we describe the low-frequency 
side of the spectra by only adjusting the mass m^(T) of 
this excitation, while keeping all other parameters of this 
mode (i.e., the high-frequency cutoff of the CM prop- 
gator, the Q coefficient, and the overall intensity coeffi¬ 
cient Ai) fixed at all temperatures, within the tempera¬ 
ture range considered here. Thus, we obtain the marked 
temperature dependence of the longitudinal NCM mass. 



Figure 5. Subtracted experimental Raman absorption spectra 
in the B 2 g channel, at various temperatures, for LSCO Oil x — 
0.04 (symbols). The theoretical fits (solid lines) consider the 
contribution of the longitudinal NCM and of the CO CM, with 
Q = 3.16 cm“^, Cc = 333 cm“^, = 50 cm“^. Ag — 7.14 

(a.u.). The other fitting parameters are reported in Fig.[^ 
The inset reports the temperature dependence of the mass of 
the longitudinal NCM (black circles) and of the CO CM (red 
squares). 


which is reported in the insets of Figs.j^ and On 
the other hand, the complete quantitative agreement be¬ 
tween data and theoretical fits is only obtainable by ad¬ 
justing more freely the secondary CO CM. This mode is 
therefore allowed to vary its parameters with T, as re¬ 
ported in Fig.[^ The temperature dependence of the CO 
CM parameters Qc and Ac oc likely reflects an increas¬ 
ing damping and a decreasing coupling to the fermion 
quasiparticles with increasing T. Of course, the estimates 
and variations of these parameters may be quantitatively 
affected if the constraint of T-independent parameters 
for the longitudinal NCM (but for its mass mg) were re¬ 
laxed. Furthermore, we cannot exclude that static ne¬ 
matic order has eventually occurred, e.g., in the sample 
with X = 0.02 at the lowest temperature. In this case our 
analysis, which is only valid above the critical tempera¬ 
ture, should be modified to deal with a broken-symmetry 
phase. This might reflect in a reduction of the primary 
peak, due to the freezing of NCM fluctuations, and could 
be the cause of the non monotonic behavior of the peak 
height as a function of T, observed in the sample with 
X = 0.02. To asses the occurrence of static nematic order 
at low temperature, a systematic experimental investiga¬ 
tion in this temperature regime is needed. 


IV. DISCUSSION AND CONCLUSIONS 

Our analysis showed that the anomalous Raman ab¬ 
sorption observed in underdoped cuprates can be inter- 
















100 




Big 

Bzg 


LSCO 

x>0.05 

CO (1,0) 

N 

N 

No nematicity at 
x>0.05: 

CO only 

LSCO 

x<0.05 

N 

CO (1,1), 

N 

Structured peak 
(shoulders): 
two CMs 

YBCO 

x<0.05 

N 

N 

No structured peak 
(no shoulders): 

N only 




Pseudogap very strong at 
x<0.05: no coupling to CMs? 


Figure 6. (a) High-frequency cutoff Qc for the CO CM for 

a sample at x — 0.02 (black empty squares) and at x — 0.04 
(red filled squares), (b) Amplitude coefficients Ac for the CO 
CM for a sample at a: = 0.02 (black empty squares) and at 
X = 0.04 (red filled squares). 


preted in terms of direct excitation of nearly critical CMs 
(see Fig.j^. The strong temperature dependence of the 
mass (i.e., inverse square correlation length) of the “pri¬ 
mary” CM, identified as the longitudinal NCM (with dy¬ 
namical critical index, z£ = 3), captures the correspond¬ 
ingly strong variation of the spectra. This CM alone fully 
accounts for the spectra of YBCO. In LSCO, instead, a 
distinct “secondary” CM, with different dynamical crit¬ 
ical index z = 2 , is needed to reproduce the composite 
lineshape. Within the two candidates considered in our 
scheme (transverse NCM and CO CM), our fits indicate 
that the CO CM is the most suitable. 

For symmetry reasons the secondary CO CM cannot 
occur in all channels: The first two rows in the sketch of 
Fig. [^summarize the findings of Ref.lbOlin LSCO as far as 
CO is concerned: The correct CO (i.e., with finite qc in 
the direction compatible with inelastic neutron scattering 
experiments) appears as an observed absorption (green 
case) only in the theoretically predicted channel. 

Two questions still remain to be answered, in order to 
complete the scheme of Fig.[^ First of all, the NCMs 
would equally contribute to the Big and B 2 g channels. 
Therefore, they would not only add to the CO fluctua¬ 
tions that give absorption in the Big channel at larger 
doping {x > 0.05) in LSCO, but would also give rise to 
absorption in the B 2 g channel. Since this is not observed 
(the corresponding box is red in Fig.[^, we infer that 
NCMs disappear in LSCO at x > 0.05 ^ee the comment 
box in the first row of Fig.[^. This is consistent with 
the observation of an increasingly stronger stripe order 
at higher doping m, where CO CM alone [along the 
(1,0) and (0,1) directions of the BZ] accounted for the 
anomalous Raman absorption at x = 0.10 and 0.12. 

The second related question is: If the NCMs are 
present and contribute to the absorption in B 2 g at low 
doping both in LSCO and YBCO, why are they not vis- 


Figure 7. Schematic comparison of the theoretical expecta¬ 
tions and the experimental observation of an anomalous Ra¬ 
man absorption. The theoretically involved CMs are indi¬ 
cated with N in the case of the NCM, while for CO we also 
report the direction of the characteristic wavevector, as estab¬ 
lished by inelastic neutron scattering. The related symbols 
only appear in the box where they are expected to contribute 
on the basis of symmetry arguments. The experimental ob¬ 
servation of an anomalous Raman absorption is depicted as a 
green case in the column of the corresponding channel. Red 
cases indicate, instead the lack of anomalous Raman absorp¬ 
tion in experiments. Our remarks and possible indications (in 
boldface) are contained in the comment boxes. 


ible in the (for them allowed) Big channel? As yet, 
we do not have a definite answer. We argue that the 
strong pseudogap occurring in lightly doped cuprates at 
T < 200 K could play a key role in suppressing the Big 
absorption. Specifically, the Big form factors select the 
quasiparticles in the fermionic loops of Fig.[2 precisely 
from the BZ regions where the pseudogap is largest. 
Therefore, only the quasiparticles in the remaining Fermi 
arcs, mostly weighted by the B 2 g form factors, remain to 
couple the Raman photons with the NCMs. In Appendix 
C, we obtained a numerical estimate of this suppression, 
finding indeed that it can be substantial. 

Based on the above discussion, we can draw the 
schematic “phase diagram” for underdoped cuprates re¬ 
ported in Fig.[^ Despite its speculative character, it 
is compatible with various theoretical and experimental 
findings, and accounts for the assessed relevance of ne¬ 
matic order in cuprates [^[38H4T] . It also complies with 
the proposal of a nematic order resulting from the melt¬ 
ing of stripes m or of a nematic or smectic phase in 
strongly underdoped LSCO (and possibly YBCO) [47], 
arising from the aggregation of doped charges in short 
segments (blue region). The orientation of these seg¬ 
ments breaks the lattice C 4 , rotational symmetry prepar¬ 
ing the route to CO at higher doping, when the segments 
merge into stripes (green region). The fluctuating char¬ 
acter of the nematic phase should give rise to nearly crit¬ 
ical fluctuations of the form of Eqs. (§ and (§. On 
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^nematicity CO^ 



Pseudogap with k-space structure ~|Bi | 

doping 

x=0.05 x=l/8 CDW-QCP 

Figure 8. Schematic evolution of the nematic (blue region) 
and stripe CO (green region) phases in underdoped cuprates 
with doping (disregarding superconductivity). The pseudo¬ 
gap region, where the Raman response in the Big channel is 
expected to be suppressed, is highlighted in red. Upon in¬ 
creasing doping, the nematic phase evolves into a CO phase, 
which in turn vanishes at a CO quantum critical point around 
optimal doping. The orientation of the segments and/or 
stripes may change with cuprate family and doping. 


the other hand, CO fluctuations become prominent by 
increasing doping and appear in the Big channel above 
X = 0.05. In LSCO these fluctuations are present (in the 
diagonal directions of the BZ) also at x < 0.05 and con¬ 
tribute to the B 2 g absorption, but the tendency of CO 
fluctuations to become more relevant at larger doping is 
clearly visible by comparing Figs.and At the same 
time, the insets of Figs.[^and|^also display an increase of 
the low-temperature limit of the correlation length of the 
CO CM upon increasing doping. Hence, nematic and CO 
fluctuations coexist in very underdoped LSCO, the pre¬ 
dominance shifting from nematic to CO with increasing 
doping. This indicates a continuous evolution from the 
nematic (charge segment) phase to the stripe phase where 
charge and spin degrees of freedom are tightly bound, 
yielding a definite relation between spin and charge in¬ 
commensurabilities (typical of the stripe phase). We re¬ 
lied on this relation to implement our symmetry-based 
selection rules for LSCO. On the other hand, our finding 
that NCMs alone are relevant in YBCO at very low dop¬ 
ing supports the idea that oriented charge segments may 
occur in this material as well, accounting for the order- 
parameter-like disappearance of the incommensurability 
in the spin response with increasing temperature mm, 
as observed in Refs. mM- The lack of CO fluctuations 
at low doping and the opposite doping dependence of the 
charge and spin characteristic wavevectors [33] indicate 
a nematic-to-CO switching different from that in LSCO. 
However, both materials seem to eventually evolve into a 
charge-density-wave phase ending into a quantum critical 
point around optimal doping, as theoretically proposed 
[4| and recently observed |33) . 


ACKNOWLEDGMENTS 

S.C. and M.G. acknowledge financial support form the 
University of Rome “Sapienza” with the Project Awards 
n. C26H13KZS9. 


Appendix A: Details on the calculation of the 
Raman response 

Inspection of Fig.[^shows that the diagrammatic struc¬ 
ture of the Raman response due to the excitation of two 
CMs involves two inequivalent fermionic loops (on the 
left of the two diagrams), which multiply two CM prop¬ 
agators and the fermionic loop on the right of the dia¬ 
grams. Calling A^^(q, cj^) the sum of the two differ¬ 
ent fermionic loops (frequencies and momenta are those 
displayed in Fig.[^, we can write the expression for the 
sum of the two diagrams as 

Xij Ay 

q,m 

^ -^A (Q? (Q? , CJ^), 

where i,j = Big,B 2 g, A ,77 = £,t,c, and Lj''{q,Pi,LOm) 
stands for the fermionic loop in the right part of the 
diagrams. The above expression can be made symmetric 
also with respect to the latter fermionic loop, observing 
that the integrated k' can be changed into — k', and the 
fermionic Matsubara frequency can be shifted to — 
ui (frequencies and momenta are those displayed in the 
fermionic loop on the right of both diagrams in Fig.[^. 
Then, exploiting the parity of the CM propagators with 
respect to both momentum and frequency arguments, we 
can take q ^ —q and ujm —^mi Summing 

the two equivalent expressions and dividing by 2 we are 
finally led to calculate 

q,m 

X Dx (q, UJjyi^Dg (q, UJ^yi + ^0Ay(q , Vi , UJjn ). 

This expression has the formal structure of a Raman re¬ 
sponse where to CMs are directly excited by the scattered 
electromagnetic radiation, and A^^(q, z//, cj^) plays the 
role of an effective Raman vertex, resulting from the sum 
of the fermionic loops with attached direct and crossed 
boson lines in Fig.[2 In Sec. m we discuss the calcula¬ 
tion of the effective vertex A-^{q^ui^ujrn)- This calcula¬ 
tion is further specialized to the cases of NCMs and CO 
CMs in Secs. mi] an d llTBll respectively. 

Appendix B: Fitting procedure of the anomalous 
Raman absorption 

To fit the anomalous contribution of the Raman ab¬ 
sorption due to two virtual CMs, as represented in Fig.[^ 




















10 


one needs to subtract the regular part of the spectra aris¬ 
ing, e.g., from the dressed quasiparticles. However, the 
subtraction procedure has to face the problem of pseudo¬ 
gap formation occurring in the underdoped regime: At 
substantially high temperatures (above 300 K) there is 
no pseudogap, which instead sets in below 200 K. The 
anomalous absorption peak we are interested in starts to 
appear on top of the (pseudogapped) broad absorption 
spectra at lower T. Thus, when the non-pseudogapped 
spectra at T ^ 300 K are subtracted from the low- 
temperature pseudogapped ones, a negative differential 
absorption is found over the frequency range of the pseu¬ 
dogap. Although this is not crucial for the qualitative 
description of the anomalous peaks, to get quantitatively 
more precise fits we exploit the fact that the pseudo¬ 
gap sets in rather rapidly and, once established, depends 
only very weakly on T. Therefore, we add a smooth 
parabolic contribution Xb — ^{^max — ^ 

with uj in cm“^ and T in K, just designed to can¬ 
cel the negative part at each temperature. For YBCO 
we take ^Imax = 1000, H(55) = 760, H(86) = 860, 
B{127) = 1100, while at T = 190, 254, 282 K no compen¬ 
sation is needed because the pseudogap is not open. For 
LSCO at X = 0.02 we take VLmax = 800, H(35) = 500, 
H(88) = 500, H(125) = 600, H(182) = 1000, while at 
T = 255 K again no compensation is needed because the 
pseudogap is not open. The same procedure is carried 
out at X = 0.04, with VLmax = 1000, B{b2) = 600, 
H(105) = 600, H(137) = 800, H(154) = 1000, H(169) = 
800, H(207) = 1500, H(252) = 2000, H(301) = 2000. 
Fig. [^exemplifies the procedure for LSCO with x = 0.04 
at T = 52 K. The blue curve represents the raw data, 
from which we subtract the red data at T = 331 K, ob¬ 
taining the purple curve with unphysical negative ab¬ 
sorption. The pseudogap effect is then corrected by the 
addition of the smooth parabolic contribution, leading to 
the final black curve. Once these differential spectra are 
thus brought to have a zero background we proceed to 
fit the strongly T-dependent anomalous peaks. 


Appendix C: Pseudogap, Fermi arcs, and Raman 
response suppression 

The strongly underdoped phase of cuprates is charac¬ 
terized by the presence of a pseudogap that strongly sup¬ 
presses the low-energy electronic degrees of freedom. In 
particular, the electronic states in the so-called antinodal 
regions of the BZ, around (Ttt, 0), (0, Ttt), are gapped, 
while the so-called nodal states, along the (±1,±1) di¬ 
rections, survive giving rise to Fermi arcs which shrink 
upon lowering temperature and doping. In this appendix 
we investigate the effects of this suppression of the low- 
energy electronic states on the coupling between the Ra¬ 
man vertices and the NCMs. Indeed the fermionic loops 
entering the diagrams of Fig.[^ involve the integration 
over fermionic degrees of freedom coupled to the nearly 
critical CMs, with the low-energy fermions being the 



Figure 9. Subtraction procedure on the raw Raman data of a 
LSCO sample at x — 0.04 and T = 52 K (blue curve and sym¬ 
bols) . The red curve and symbols correspond to the raw data 
at T = 331 K. Once the latter are subtracted from the former, 
the purple curve and symbols are obtained. To eliminate the 
negative absorption, the parabola Xh — ^(1000 — Ct;)/500^, 
with UJ in cm“^, is added to finally yield the absorption re¬ 
ported with the black curve and symbols. 


most effective in coupling to the low-energy CMs. There¬ 
fore the opening of gaps in the electronic spectra nat¬ 
urally entails a substantial reduction of the overall re¬ 
sponse of the CM. However, the Raman vertices 7 i(k) 
weight differently the fermionic states along the Fermi 
surface and it is therefore quite natural that the fermionic 
loops are differently suppressed depending on the chan¬ 
nel. To estimate this effect is the aim of this appendix. 
More specifically, we will consider the NCMs only, be¬ 
cause the CO modes in the very underdoped LSCO were 
shown in Ref. na to be only visible in the B 2 g channel. 
So it would be meaningless to compare the pseudogap 
effects in the two Raman channels. On the contrary, the 
NCM are singular at q ~ 0 and therefore should give a 
strong contribution to the Raman response in both chan¬ 
nels. This appendix will instead demonstrate that the 
interplay of Raman vertices and momentum dependence 
of the pseudogap strongly suppress the loop in the Big 
channel in comparison to the B 2 g case. 

We adopt the simplifying assumption of a circular 
Fermi surface and we start from the expression for the 
fermionic loop, Eq. (H). Since the NCMs are singu¬ 
lar at small q, we expand in this limit the quantity 
— ^k+q- Expanding up to order the de¬ 
nominator and exploiting the S function to perform the 
integral along the radial momentum variable, one obtains 


A, 


2M^ 


p27r 

/ deF{e)- 

Jo 1 


7i(6') 


_q^ 

2kl 


■ cos(26> — 20) 


(Cl) 


where 0 is the angle between q and the x axis in recip¬ 
rocal space. At this stage we have phenomenologically 
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Figure 10. Square of the fermionic lo ops as a function of the 
angle (j) calculated according to Eq. ( |C1[ ) in the Big channel 
(solid lines) a nd i n the B 2 g channel (dashed lines). The F{0) 
function [Eq.( |C2| ] has been set to produce arcs shrinking the 
Eermi surface by a factor 1 (©m = 7r/4, whole Eermi surface, 
black curves), 0.75 (©m = Stt/IG, red curves), 0.5 (©m = 
tt/S, green curves), and 0.25 (©m = 'tt/IG, blue curves). The 
parameter = 0.05, was set to smoothen the angular cutoff 
producing the arcs. The boson momentum q has been chosen 
such that I2k\ — 0.01, while M — 


introduced a function 

H , {[0-(2n-l)7r/4]2-e2^}/A= ’ 

n 

with n — 1, 2, 3,4, which simulates the effect of the pseu¬ 
dogap on the Fermi surface of the 1 — 4 quadrants. Specif¬ 
ically this function leaves the states near the diagonal un¬ 
touched, while for ©m < 7r/4 it rather sharply suppresses 
the integration in the gapped antinodal regions for ^’s 
far from the nodal direction ^ = 7r/4 (for ©m = 7r/4 one 
recovers the full ungapped Fermi sur face) . This essen¬ 
tially restricts the integration in Eq. ( |Cl| to the angles 
of a Fermi arc allowing to explore the different action of 
the Fermi surface shrinking on the value of the fermionic 
loop. The parameter measures how rapidly the pseu¬ 
dogap is switched on and off along the FS, and we take 
it to be much smaller than 7r/4. 

Fig. [^displays the square of the fermionic loops in the 
two Raman channels as a function of the angle (\) between 
the boson transferred momentum q and the x axis. The 
calculation clearly shows the increasingly strong suppres¬ 
sion of the B\g fermionic loop (solid curves) upon reduc¬ 
ing the length of the Fermi arcs. The suppression is much 
less pronounced in the B^g fermionic loop (dashed lines). 
These results are rather natural because the pseudogap 
suppresses the states that more effectively contribute to 
the Big loop, while the fermion quasiparticles contribut¬ 
ing more to the B 2 g channel survive in the Fermi arcs. 
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